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ABSTRACT 



The existence of millisecond pulsars with planet-mass companions 



( iBailes et al. 



20111 ) in close orbits is challenging from the stellar evolution point 
of view. We calculate in detail the evolution of binary systems self-consistently, 
including mass transfer, evaporation and irradiation of the donor by X-rays feed- 
back, demonstrating the existence of a new evolutionary path leading to short 
periods and compact donors as required by the observations of PSR J1719-1438. 
We also point out the alternative of an exotic nature of the companion planet- 
mass star. 



Subject headings: pulsars: general — pulsars: individual (PSR J1719-1438) — stars: 
evolution — binaries: close 
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INTRODUCTION 



The recent report (IBailes et al 



20 111 ) of a binary millisecond pulsar in a 2.2 h orbit 
featuring a Jupiter-like mass companion with a lower bound for the mean density of 
p > 23 g cm~ 3 is both important and challenging for stellar evolution theory. Indeed, 
the role of the pu lsar wind and illumination feedback have been deemed as important 



(IBailes et al. 



2011 



), but it was not clear whether the interplay of all the effects is enough to 
reproduce the observed features. In order to account for the existence of the PSR J1719-1438 
system we have looked for an evolutionary scenario in which a normal star evolves losing 
most of its mass and reaching the observed configuration. We have considered close binary 
systems composed by an accreting neutron star (NS) orbiting together with a normal donor 
star. We attempt to answer the full history of these systems below, and report the first 
complete results in this work. 



2. CALCULATIONS AND RESULTS 



In the case of this kind of systems, the donor star evolves essentially as an isolated 
object up to the moment at which its radius R2 nearly equals the radius of the Roche Lobe 
Rqj. This phenomenon is usually referred to as the onset of the Roche Lobe Overflow 
(RLOF). Near RLOF, tidal dissipation forces the orbit to become circular with an period 
Pi, a starting point for the calculations. For the case of low mass donor stars with masses 
fulfilling the condition < M2/M1 < 0.8, the radius of a sphere with the volume of the 



Roche Lobe can be approximated by (jPaczyhski 



197lh 



l As usual, we shall refer the NS (donor star) as the primary (secondary) with subindex 
1(2). 
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Rt 



0.46224 a 



Mo 



M 1 + M 2 



1/3 



where a is the semiaxis of the circular orbit. Since that moment on, the donor star 
transfers mass across the Lagrangian point L\ towards the NS. This process, in turn, 
makes the orbit to evolve. At present it is not clear how much of the matter transferred 
from the donor star is effectively accreted by the NS. Hereafter we define /3 as the fraction 
of transferred material that is accreted by the NS ( M-\ = —(5Mi, but always b elow the 



Podsiadlowski et al. 



20021 ). If < 1 



Eddington limit M Edd = 2 x 10~ 8 M Q yr' 1 ; see, e.g., 
some material is lost by the system carrying away the specific angular momentum of the 
secondar y. As the value of is not critical in determining the evolution of this kind of 



systems (IDe Vito & Benvenuto 



and more recently in p opulation synthesis calcu lations by 



gravi tational radiation (ILandau fc Lifshitz 



i Podsiadlowski et al. 


f 


2002), 


Belczvnski et al. 


( 


2008 


). Also, 


netic braking (IVerbunt & Zwaan 



198fl ) are known to provide relevant angular momentum sink 



lie LIU I 

inkfl 



M 2 due to RLOF is described by the expression given by 



Ritterl (|1988h . 



M 2 .h>L(>r --Uioxp ( 2 H Rl 



(2) 



where M is a smooth function of M-\ and M ?, and Hp is the pressure scale height at the 



photosphere (for further details, see 



Ritter 



19881 ) . The above given description corresponds 



to t he standard treatmen t for t he evolution of Low Mass X-Ray Binaries (LMXBs); see, 



e.g.; 



Podsiadlowski et al. 



( 120021 ). However, as it will be clear below, these ingredients are 



2 Another law of magnetic braking has been presented in llvanova fc Taaml (120031 ); adopt- 
ing it the results of the present work may change. This will be explored in a future paper. 
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not enough to account for the formation of a binary pair like PSR J1719-1438. Another 



effect that drives further mass 
the pulsar radiation. Following 



oss from the donor star is the evaporating wind, driven by 



Stevens et al. 



( 119921 ). we include this effect by considering 



M< 



2,evap 



2v 2 



V a 



2, esc 



(3) 



where the pulsar's spin down luminosity Lp is given by Lp = 47r 2 J 1 P 1 P 1 (Ji is the 
moment of inertia of the NS, Pi is its spin period and Pi its spindown rate), v^^sc is the 
escape v elocity from the dono r star surface and / is an efficiency factor that will be set to 



0.1 as in 



Stevens et al 



(|1992| ). As we shall be concerned with very short orbital periods, 
a relevant phenomenon to be considered is irradiation feedback. When the donor star 
transfers mass onto the NS, it releases an accretion luminosity that illuminates the donor 
star with a flux P,> r = t^^h^Mi, where a,i rr is a constant that accounts for the fact that 



i-KO 2 Rl 

neither all the lumino sity has to be re 



emitted isotropically terming fc Ritter 



eased as electromagnetic radiation, nor has to be 



20041 ) . The radiation incident onto the donor star 



partially blocks the rel ease of its internal 



has been addressed by 



energy, modifying its evolution. This problem 



Hameury fc Ritterl (I1997I ). Here we shall assume the validity of the 



point-source model, i.e., that the accreting NS is the only source of radiation incident onto 
the donor star. 

In order to compute the evolution of these systems we have emp l oyed o ur 



Benvenuto fc De Vitol ( 120031 ) and 



detailed (Henyey) evo l ution ary code described in 

De Vito Sz Benvenuto I2OI2 ). which has been modified to incorporate irradiation feedback 
and evaporating winds as described above. For a compact binary system to be an adequate 
candidate to account for the properties of PSR J1719-1438, it must have a very close orbit. 



The RLOF will occur during the h ydrogen core burning stage; t 



as a Class A mass transfer episode (jKippenhahn &: Weigert 



lis is usually classified 



19671 ) (whereas Class B and 
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C episodes correspond to the cases in which the onset of the RLOF occurs after the 
exhaustion of core hydrogen and helium respectively). An exploration of the parameter 
space defining a particular compact binary system (M 1; M 2 , P») indicates that there exists 
a restricted region of which that leads to the formation of systems like PSR J1719-1438. 
For example, there is a narrow range for the initial orbital period: if P, is too short (say, 
< 0.5 d) even at the minimum radius (on the ZAMS), the donor star would be transferring 
mass. On the other side, if Pj is larger than about « 0.9 d, the system evolves on an orbit 



that widens enough to allow for the formation of a low-mass (~ 0.25M^) he 



dwarf (He WD) - millisecond pulsar pair (see, e.g., 



De Vito & Benvenuto 



ium white 



20121). Thus, the 



values of Pj leading to objects on converging orbits is very restricted. If the formation of 
the PSR J1719-1438 system proceeded the way we considered here, P, should have fallen 
in this interval^. Furthermore, as we shall see below, the system has to evolve for a quite 
a long time to reach mass values as low as those indicated by observations. This, in turn, 
imposes a lower limit for its initial mass: if it is very low, the system would need to evolve 
for a time in excess of the age of the Universe. For more massive stars, there exists a limit 
im posed by the stability of th e mass transfer at the onset of the RLOF (for further details 



sec 



Podsiadlowski et al. 



20021 ) . This set of conditions strongly suggests that these systems 



should be rare. 



Observations of the pulsar signal are quite s table, and thus d o not suggest that 



PSR J1719-1438 is undergoing a RLOF episode ( IBailes et al. 



2011 



M. Bailes, private 



communication). Thus, the donor star should be smaller than its corresponding Roche 



3 In any case, here we should remark that the interval of Pj referred above depends on 
the particular physical ingredients assumed in our computations that are certainly not fully 
known. The related present uncertainties should affect the precise value of the period interval, 
although it will be still within a narrow range. 
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Lobe. If RLOF were the only process giving rise to mass loss from the donor star, this 
fact would be very difficult to account for. It is well known that low mass WDs behave 
like polytropic spheres of index n = 1.5. For these structures, the mass-radius relation 
is R oc M -1 / 3 . Therefore, the star expands in response to mass loss. This behavior 
continues as long as the equation of state is dominated by electron degeneracy. As the 
star experiences further mass loss, it lowers its density and the degree of degeneracy; and 
non-ideal effects become more important. Eventually, the mass-radius relation changes, and 
there is a mass value for which the radius passes through a maximum. For example, if the 



object has a helium dominated composition this corresponds to an ob 
M?, = 2 x 10 _3 M o , which has a radius of R 2 = 5 x 1O~ 2 -R (see, e.g., 



ect with a mass o 



Deloye fc Bildsten 



20031 ). As the less massive object in the system losses mass (and angular momentum), we 
arrive to a situation in which the orbit gets wider while keeping R 2 — Rl ~ Hp. When the 
donor star reaches the maximum-radius mass value, a further mass loss will force the star 
to contract, detaching the donor star from its Roche lobe (notice that for these mass values 
of the donor star, the timescale of orbital evolution due to gravitational radiation is too 
long to lead the donor star into contact). 

If mass loss/transfer were only due to RLOF episode(s), in reaching the observed 
configuration, the system would need a timescale in excess of the age of the Universe. A 
natural way out of this apparent paradox is provided by the evaporating wind described 
above. As a matter of fact, during the advanced stages of evolution in which the donor star 
mass becomes very low (M 2 < 2 x 1O _2 M ), such evaporating wind dominates the donor 
star mass losses and the orbital evolution, even if the system is still on RLOF conditions. 
This effect makes the orbit to become wider than it would be if we consider RLOF solely, 
making the star to detach from its Roche lobe before reaching mass values as low as 
that corresponding to the maximum radius for its composition. In order to explore the 
plausibility of this scenario we have computed the evolution of several systems assuming 
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a solar composition donor star with an initial mass value of M2 = 2M , a "canonical" 
NS of Mi = 1.4M , and some values for the initial orbital period that lead to this kind 
of binary systems Pj = 0.75 d, 0.80 d and 0.85 d. We considered evolutionary sequences 
with and without irradiation feedback. In this paper we shall not discuss the process of 
formation of main sequence-NS CBSs from which we star out our calculations. Also, we 
should warn the reader that it is expected to be possible to arrive to a configuration like 
that of PSR J1719-1438 from other initi al conditions, different from the ones we assumed. 



These proces ses have been discu s sed b y 
The work by 



van Haaften et al. 



Belczynski fc Taaml (120041 ) and references therein 



( 120121 ) noticed several problems in the formation of the 



system and attempted to model the outcome varying the donor and its wind. 

In Fig. [TJwe show the orbital period of the system as a function of the donor mass. It is 
remarkable that irradiation feedback does not induce any dramatic effect on such a relation. 
The observed period for PSR J1719-1438 indicates that for each model two solutions exist: 
one with M2 > 0.10M Q , (while the orbit is shrinking) and other with Mi = O.OIM^ (w hile 



the orbit is expanding) . In view of the mass function of this system ( iBailes et al. 



201 If ) 



(where P is the orbital period, is the semiaxis of the pulsar orbit and i is the 
inclination angle of the orbit with respect to the line of sight), it is difficult to consider the 
first solution as physically plausible. For M\ = 1.4OM , we would need sini w 0.01, which 
has a very low probability, whereas for the other solution we still need small but tolerable 
values of sini w 0.1. Notice that because at late times the evaporating wind dominates 
mass loss/transfer (see below), it accelerates the evolution of the system (as compared with 
the standard case in which this effect is ignored) making it possible to reach a configuration 
compatible with the observed state of PSR J1719-1438 within a long but acceptable 
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timescale of 6-7 Gyr. We show in Fig. [2] the mass transfer rate for the case of Pi = 0.8 d, 
with and without irradiation feedback. Due to irradiatio n, the donor star u nderg oes cyclic 



mass transfer episodes in a way similar to that found by 



Biining fc Ritterl ( 120041 ) . Notice 



that this oscillating behavior is restricted to an intermediate stage of evoluti on. Despite t he 



2008|), 



uncertainties associated with the present treatment of irradiation feedback (iRitterl 
it is a fortunate situation that the final properties of the system are largely independent 
of the former. We show in Fig. [3] the evolution of the donor and the Roche lobe radii, 
demonstrating that the system ultimately detaches around 6 Gyr. Finally, in Fig. [4] we 
show the evolution of the mean density of the donor star p. We find that since ages ~4 Gyr 
on (well before detachment from the Roche lobe), p overcomes the lower limit deduced from 
observations. 

Regarding the final internal composition of the donor star, it depends on the value of 
Pi. For the shortest possible initial orbital periods, core hydrogen burning is quenched by 
mass loss (internal temperature falls down fast enough to appreciably slow down nuclear 
activity) and the final hydrogen abundance is 0.45 by mass. For the largest Pi for which 
CBSs evolve to a black widow configurations, compatible with the characteristics of the 
PSR J1719-1438 system, hydrogen is almost absent being a helium-dominated composition. 
It is worth to remark that for the formation path we addressed in this paper, no other 
composition is possible for the donor star interior. 



3. DISCUSSION AND CONCLUSIONS 

The conclusion of this study is to identify a definite new path for the evolution of 
binary systems evolving into planet-like - millisecond pulsars pairs, featuring R 2 < Rl 
for the donor star and a mean density p > 23 g cm~ 3 for it. Our calculations show 
self-consistently that this is indeed possible, even for objects composed by a mixture of 



- 10 - 



hydrogen and helium, without the need of postulating a carbon interior, on a reasonable 
timescale. The initial conditions for this evolution are actually quite stringent, as identified 
above; otherwise the outcome of the evolution is very different. 

Finally, the exciting possibility that the companion of the millisecond pulsar 
PSR J1719-1438 is not a WD-like but a truly exotic object (i.e. composed of some form 
of quark matter) should not be overlooked. This would easily explain why there is no 
modulation even for an edge- on inclinatio n . Act ually, in the latter case the stringent 



photometric limits derived in iBailes et al. 



( 1201 ll ) using Keck-LRIS instrument cannot 



be used to place constraints to the inclination, because the absence of signal is a quite 
natural outcome. In other words, for the cases of a strange quark matter nugget or 
structured strangelet chunk, the size of the companion would be too small to detect any 
photometric signal. The exotic model also predicts that no carbon/helium lines should be 
ever observ ed associated to th e companion. In addition, the lack of detection of evaporation 



signatures (IBailes et al. 



201 ll ) would be naturally accommodated. Finally, and because 



of angular momentum considerations, we expect that the orbit angular momentum J of 
the quark compan i on to be aligned with the spin of the pulsar "born in original spin", 



e.g., 



Camilo et al. 



( 11994] ) . These are quite strong, albeit straightforward predictions to be 



checked in futures studies addr essing the nature o f this system. The p roposals of extended 



exoti c stars (strangelet dwarfs, 



Alford et al. 



2012 



and strange dwarfs, 



Glendenning et al. 



19951 ). would need an evaluation of their surface properties, which would still depend on 
the existence or absence of a normal matter atmosphere to reprocess the incident pulsar 
radiation. This would be difficult to distinguish from conventional helium or carbon WDs. 
However, in these scenarios there is no link between the evolution of the system and the 
final masses and period, and the millisecond pulsar could be very young and not recycled 



at all. 
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Fig. 1. — The orbital period-mass relation for the donor star corresponding to systems with 
a solar composition, 2M & normal star and a 1.4M Q neutron star in orbits with initial periods 
Pi of 0.75 d (black), 0.80 d (red), and 0.85 d (blue) respectively. Full lines correspond to 
calculations neglecting irradiation feedback, while the results for 0.80 d and two values for the 
irradiation parameter «j rr (0.01 -cyan- and 0.05 -pink-) are shown with dashed lines. Stars 
spend 1 Gyr evolving leftwards from one mark to the next one along the trajectories. The 
observed orbital period for PSR J1719-1438 is marked with an horizontal dashed line. These 
systems attain the observed period with adequate masses (< O.O5M ) after long (6-7 Gyr), 
but acceptable timescales. 
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Fig. 2. — The evolution of the mass transfer rate from the donor star for the case Pj = 0.8 d, 
ignoring (black line) and considering (blue line) irradiation feedback (with ai rr = 0.05) as in 
Fig. [TJ The mean value of the mass transfer rate is very similar, irrespective of the inclusion 
of irradiation feedback. However, for ages between 2 and 4 Gyr, irradiated mode ls undergo 



Buning fc Ritterl fl2004h . Some of 



a sequence of RLOF episodes similar to those found by 
these episodes are depicted in further detail in the inset A where black, red and blue lines 
represent the total mass transfer rate, the RLOF and evaporation contributions respectively. 
Finally, in inset B we show the same quantities (with colors having the same meaning as in 
inset A) at the end of the non-irradiated sequence (irradiated ones behave in a similar way). 
Remarkably, after 5.5 Gyr mass transfer rate is dominated by the evaporation wind driven 
by pulsar irradiation. 
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Fig. 3. — The evolution of the radius of the donor star R 2 (solid lines) and its corresponding 
Roche lobe Rl (dashed lines) for the same models of Fig. [2, ignoring (black lines) and 
considering (blue lines) irradiation feedback (with a irr . = 0.05). In the case of the irradiated 
sequence, the donor star suffers from a series of contractions and expansions corresponding 
to the cyclic mass transfer regime, shown in detail in the inset. However, the long term 
evolution of the radii of irradiated models is very similar to those of non-irradiated ones. 
Points labelled with A and C correspond to the stages at which the orbital period is equal 
to the observed value, while B indicates the final detachment of the donor star. Notice that, 
remarkably, due to the evaporation wind driven by pulsar irradiation, the donor star attains 
the orbital period observed for PSR J1719-1438 with low mass values (M 2 ~ 0.010M Q ) in 
detached conditions, as indicated by observations. 
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Fig. 4. — The evolution of the mean density for the donor star for the same sequences shown 
in Fig. [TJ Lines and colors have the s ame meaning as the re. We show with a horizontal line 



the minimum mean density found in 



Bailes et al. 



( 120111 ) inferred from observations. First 



the systems attain the observed orbital period when under RLOF conditions, when the mean 
density is lower than the value p = 23 g cm~ 3 . After that the star detaches from the Roche 
lobe and the mean density increases above the referred minimum value, in agreement with 
observations. 
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Table 1. Parameters of the P SR J1719-1438 emp loyed in the calculations (from 



Bailes et al. 



(1201 lh 



Parameter 



Value 



Epoch (MJD) 

Port (d) 

a p sini (It-s) 
p(g cm -3 ) (inferred) 



172.70704459860(3) Hz 
-2.2(2) X lO" 16 

55411.0 
0.090706293(2) 
0.001819(1) 
> 23 



